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ABSTRACT 

The mechanisms of formation of a metastable defect isomer of fullerene Cgo due to the 
Stone- Wales transformation are theoretically studied. It is demonstrated that the paths 
of the "dynamic" Stone- Wales transformation at a high (sufficient for overcoming potential 
barriers) temperature can differ from the two "adiabatic" transformation paths discussed 
in the literature. This behavior is due to the presence of a great near-flat segment of the 
potential-energy surface in the neighborhood of metastable states. Besides, the sequence of 
rupture and formation of interatomic bonds is other than that in the case of the adiabatic 
transformation. 
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Various models of fullerene growth have been proposed [1]. Because a sphere-shaped clus- 
ter C 60 is formed at a high temperature under actual operating conditions, its structure can 
be significantly different from the structure of ideal buckminsterfullerene with 1^ icosahe- 
dral symmetry. On the "surface" of buckminsterfullerene, carbon atoms are arranged at the 
vertexes of 20 hexagons and 12 pentagons, which are isolated from each other. Therefore, 
regardless of the mechanism of formation of such a spheroidal cluster from graphite fragments 
and/or carbon dimers, the question arises of the paths of its evolution into an equilibrium con- 
figuration of buckminsterfullerene, that is, the question of the mechanisms of defect annealing 
that reduces the potential energy E pot of the cluster. 

Defect annealing occurs through local rearrangements of C-C bonds and is the reverse 
process for defect formation. The so called Stone- Wales transformation [2] , which consists in 
the rearrangement of two C-C bonds in buckminsterfullerene (see Fig. 1), results in a defect 
isomer closest to buckminsterfullerene in terms of energy. This metastable isomer exhibits C 2v 
symmetry and contains two pairs of pentagons with common sides. Among a great number 
of other isomers, this isomer stands out as the last segment (before buckminsterfullerene) in 
the chain of sequentially decreasing energies of a Ceo cluster on defect annealing [3]. 

In this context, attention was focused on studying changes in the mutual arrangement of 
atoms in the Stone- Wales transformation and on determining the height U of a minimum 
energy barrier to this transformation [4-10]. Two different transformation paths were con- 
sidered: (I) the rotation of a C-C bond shared by two hexagons through 90° so that all of 
the atoms remained on the cluster "surface" and (II) the rupture of a C-C bond shared by 
pentagon and hexagon; thereafter, one of the atoms initially rose above the cluster " surface" 
and then relaxed to form new bonds. In accordance with the calculations of E pot as a function 
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of atomic coordinates {Ri} performed using the density functional theory [5 - 10], the barrier 
height for path I was C// = 6-8 eV, whereas U u = 6.64 - 7.6 eV for path II. A considerable 
scatter in the values of Uj and Un does not allow one to determine unambiguously which 
of the two paths exhibits a lower barrier and, consequently, by which mechanism the Stone- 
Wales transformation occurs. Even with the use of the same exchange-correlation functional 
and equal basis-function sets (or equal cutoff energies in a plane-wave basis) for calculating 
Ui and Un, different results were obtained: Uj = Un = 7.27 eV [10], Ui = 6.30 eV < Un = 
6.64 eV [9], and Uj = 8.1 eV > U n = 7.6 eV [6, 7]. 

Note that the determination of the height of a barrier and the type of a transition state 
based on the analysis of the shape of a potential-energy surface implies an adiabatic transition 
of the cluster from one state to another. That is, in such a transition, one or more properly 
chosen interatomic distances change in accordance with a certain law with a small step along 
the reaction coordinate, whereas relaxation to a minimum (with consideration for imposed 
restrictions) total energy is performed at each step in all of the other degrees of freedom. 
Physically, an adiabatic transition corresponds to a forced deformation of the cluster at T = 0, 
and it must not occur identically to transitions at a finite (especially, high) temperature. 
In this work, we studied Stone- Wales transformation paths both in an adiabatic transition 
and at a high temperature (sufficient for overcoming a potential barrier) using a molecular 
dynamic simulation of the "life" of buckminsterfullerene. A comparison between the results 
allowed us, on the one hand, to find common features in the low-temperature and high 
temperature buckminsterfullerene metastable isomer transitions and, on the other hand, to 
reveal considerable differences between them. 

We used a tight binding potential for calculations. This potential incorporates the contri- 
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bution of all the valence electrons (four electrons from each atom) to the total energy in an 
explicit form, and it is actually an n— body potential, where n is the number of atoms in the 
cluster. Because of this, this potential surpasses simple but less reliable empirical potentials. 
Although this method is not as strict as ab initio approaches, it is highly competitive with 
them in the accuracy of the description of carbon systems [11]. Moreover, this method con- 
siderably facilitates the simulation of dynamic processes because it is not resource-intensive. 
Previously, we used this method for an analysis of the thermal stability of metastable Cs and 
C 2 o clusters [12 - 14]. For buckminsterfullerene, it gave the binding energy E b = 60£'(Ci) - 
E(Cq ) = 6.86 eV/atom and the bond lengths I = 1.396 and 1.458 A, which are consistent 
with the experimental values of = 6.97 - 7.01 eV/atom [15] and I = 1.402 and 1.462 A 
[16]. The HOMO-LUMO gap A = 1.62 eV is also consistent with the experimental value of 
A = 1.6 - 1.8 eV [17]. 

Initially, we studied the mechanisms of the buckminsterfullerene — > metastable isomer 
isomerization by adiabatic transitions. For this purpose, we calculated the potential relief 
E pot ({Ri}) for the Ceo cluster in the vicinity of an equilibrium atomic configuration (buck- 
minsterfullerene) and the closest in energy metastable atomic configuration (metastable iso- 
mer) for the above two Stone- Wales transformation paths. Figure 2 shows the results of this 
calculation. The energy of the metastable isomer is higher than the energy of buckminster- 
fullerene by AE = 1.42 eV, which is consistent with the value of AE = 1.4 - 1.7 eV obtained 
by the density functional method [5 - 10]. In both path I and path II, the dependence of 
E pot on reaction coordinate X exhibits the only stationary point (a maximum), which is a 
saddle point for the E pot ({R i }) surface and determines the barrier height. According to our 
calculations, the value of Uj = 6.46 eV is somewhat lower than U H = 6.58 eV. 
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Figure 3a shows the atomic configuration of the transition state for path I. It is symmetric 
in accordance with available published data [6 - 10]. On the contrary, the transition state 
for path II is asymmetric (Fig. 3b). One of the atoms is arranged over the cluster "surface", 
and it forms two bonds (sp hybridization). In this case, another ("surface") atom has four 
neighbors (sp 3 hybridization). This transition state has been reported previously [4, 6, 7, 9, 
10]. Note that path II exhibits an almost horizontal segment near a maximum in the E pot (X) 
curve. The atomic configuration that corresponds to the value of X at which the second 
derivative d 2 E pot (X)/dX 2 exhibits a minimum (Fig. 2, point 4) is visually similar to the 
transition state (Fig. 3b). Murry et al. [6, 7] stated that this configuration corresponds to a 
local minimum of E pot (that is, actually, to another metastable state). In recent publications 
[9, 10], no additional minimum of E pot was reported. 

Nevertheless, we found that both of the above Stone- Wales transformation paths are char- 
acterized by barriers with approximately equal heights: Uj — 6.46 eV and U n = 6.58 eV. 
These data are consistent with first-principle calculations [5 - 10], which gave similar values 
of Ui and Un in a range of 7 ± 1 eV. One might expect that the Stone- Wales transformation 
under conditions of a real experiment (that is, due to thermally activated processes) occurs 
via both of these paths with approximately equal probabilities. To test this hypothesis, we 
numerically simulated the dynamics of buckminsterfullerene at T = 4000 - 5000 K and a fixed 
total energy of the cluster, E = E kin + E pot = const. The cluster temperature was found from 
the equation [18] k B T(3n — 6) = E kin , where n = 60, and E kin is the kinetic energy in the 
system of the center of gravity averaged over 10 3 steps of molecular dynamics (the single step 
time t = 2.72 ■ 10 16 s was about 0.01 of the period of the highest frequency mode of cluster 
vibrations). On this formulation of the problem, the temperature T is a measure of energy of 
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relative atomic motion [18]. In this case, the value of E kin = 30 - 40 eV is sufficient for over- 
coming the potential barriers Ujjj ~ 6.5 eV (although E kin » Ujjj, buckminsterfullerene 
isomerization processes at T pa 4500 K occur rarely, approximately once a time of 10 6 to)- 

We analyzed in detail 45 buckminsterfullerene — > metastable isomer transitions and vice 
versa. We found that the Stone- Wales transformation occurred via path I (or by a visu- 
ally very similar mechanism) in the majority of cases (37), see Fig. 4; it occurred via 
path II only in three instances, see Fig. 5. Thus, although both of the paths are almost 
equally favorable in terms of energy (see Fig. 2), path I is much more preferable from a 
dynamic standpoint. We believe that this is due to a dramatic difference between the fre- 
quency factors Aj and An in the Arrhenius temperature dependence of the rates of transition 
kjjj = Ajjjexp(Uijj/k B T). Consequently, in the consideration of Stone-Wales transfor- 
mation mechanisms, attention should be focused on an analysis of the frequency factors of 
transitions via these two paths rather than on a further increase in the accuracy of the cal- 
culation of barrier heights Ujjj (and on determining which of these barriers is lower even by 
fractions of eV). 

In addition to the transformation of buckminsterfullerene via paths I and II, we also ob- 
served transitions accompanied by the formation of either a large "window" or three adjacent 
windows and the appearance of a branch of three atoms, which was perpendicular to the 
cluster "surface". All of these intermediate configurations correspond to potential-energy 
surface segments above the saddle points of transitions via paths I and II, that is, to other 
transformation paths with higher barriers. At T = 5000-5500 K, the relative number of these 
transitions increases dramatically. They can occasionally result in the formation of isomers 
with higher energies (for example, fullerenes with three or more pairs of adjacent pentagons 
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or nonclassical fullerenes with tetragons and heptagons). 

Note that, in the dynamic transformation of buckminsterfullerene via path I, as a rule, 
one or two atoms in the C-C bond that rotated through 90° (Figs. 1 and 3a) rose above 
the cluster "surface" through different heights every time as in the transformation via path 
II (but, unlike path II, without the formation of an sp 3 configuration). This suggests that 
there are a great number of configurations with energies close to the saddle-point energies of 
paths I and II. That is, in other words, a segment in the potential-energy surface between two 
saddle points is almost planar. To test this hypothesis, we calculated E pot as a function of 
two reaction coordinates H and and found that, indeed, the E pot (H, 0) function exhibits 
a very great radius of curvature along the coordinate H (Fig. 6). Note that a second-order 
stationary point (characterized by two imaginary frequencies) occurred between the saddle 
points of paths I and II. The energy of this stationary point is 0.56 eV higher than the energy 
of the saddle point of path II. This stationary point corresponds to a configuration in which 
one of the atoms of the rotated C-C bond rose over the cluster "surface" (however, to a lesser 
degree than in the transformation via path II). 

In conclusion, let us consider the rupture and formation of interatomic bonds in the dynamic 
Stone- Wales transformation. Figure 4 shows a sequence of snapshots of the cluster in the 
course of the transformation via path I. It can be seen that the transformation began with the 
rupture of a single bond rather than the concerted rupture of two bonds shared by pentagons 
and hexagons, as in an adiabatic case [10]. Physically, this is easy to understand: a distortion 
of the ideal buckminsterfullerene structure due to the thermal vibrations of constituent atoms 
caused a symmetry breakdown in the arrangement of these two bonds. Because of this, one of 
these bonds became weaker than the other, and it was ruptured initially. Subsequently, the 
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second bond was ruptured. Thereafter, a configuration similar to the transition state in an 
adiabatic case was formed (Fig. 3a). At the next step, new bonds were formed (consecutively 
rather than simultaneously) to complete the transformation. Note that the overall process 
occurred very rapidly and took only At/ 86 fs. The transformation via path II began 
with the rupture of either a bond shared by pentagon and hexagon or a bond shared by two 
hexagons. For the latter case, Fig. 5 shows a sequence of snapshots of the C 60 cluster. The 
configuration corresponding to the transition state (Fig. 3b) can be clearly seen in Fig. 5. 
The transformation time At// 354 fs is much longer than the value of At/; this is consistent 
with a greater "length" of path II (see Fig. 6). 
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Fig. 1. Formation of two pairs of pentagons with common sides on the rearrangement 
of two C-C bonds in fullerene Ceo (Stone- Wales transformation). For clarity, atoms in the 
background are not shown. 
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Fig. 2. Potential energy E pot of the C 60 cluster vs. the reaction coordinate X in the 
neighborhoods of (1) equilibrium and (2) metastable atomic configurations in the Stone- 
Wales transformation (Fig. 1). The energy of 60 isolated carbon atoms was taken as zero 
energy. The solid and dashed lines indicate paths I and II, respectively (see the text); points 3 
are E pot (X) maximum points (saddle points for E pot ({Ri}), and ^ is a d 2 E pot /dX 2 minimum 
point for path II. The path length along a trajectory that passes through the corresponding 
saddle point in (3n6)-dimensional space and joins buckminsterfullerene with the metastable 
isomer was chosen as the reaction coordinate. X = at saddle points. 



12 



(a) (b) 




Fig. 3. Atomic configurations of transition states 3 in Fig. 2 in the Stone- Wales trans- 
formation via (a) path I (symmetric transition state) and (b) path II (asymmetric transition 
state). Atoms in the background and C-C bonds more than 2 A in length are not shown. 
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Fig. 4. Snapshots of the C 6 o cluster in the course of the dynamic Stone- Wales transforma- 
tion via path I. Initial temperature Tj„j = 4465 ± 5 K. The first configuration corresponds to 
the time t' = 0.245 ns, and the subsequent configurations correspond to the times t' + At, 
where At = 5.4, 17.7, 69.9, 73.4, and 85.7 fs, respectively. Atoms in the background and C-C 
bonds more than 2 A in length are not shown. 
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Fig. 6. Potential energy E pot of the C 60 cluster vs. two reaction coordinates H and 6 in 
the neighborhoods of (1) equilibrium and (2) metastable atomic configurations. The energy 
of 60 isolated carbon atoms was taken as zero energy. The solid and dashed lines indicate 
paths I and II, respectively; 3 indicates corresponding saddle points; and S is a second-order 
stationary point (see the text). The reaction coordinates were chosen as follows: H is the 
height through which an atom of the rotated C-C bond raised over the cluster "surface" 
normalized to a maximum height through which this atom rose in its motion along path II 
(H = for path I; a maximum value of H — 1 was at the saddle point of path II) and is the 
angle of rotation of the C-C bond normalized to 90° (0 = for buckminsterfullerene, = 1 
for the metastable isomer, and 0.5 for the saddle points and the second-order stationary 
point S). 
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